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We present a detailed investigation of the momentum-dependent self-energy at zero fre- 

quency of weakly interacting bosons at the critical temperature T c of Bose-Einstein condensation in 
dimensions 3 < D < 4. Applying the functional renormalization group, we calculate the universal 
scaling function for the self-energy at zero frequency but at all wave vectors within an approximation 
which truncates the flow equations of the irreducible vertices at the four-point level. The self-energy 
interpolates between the critical regime k <C k c and the short-wavelength regime k 3> k c , where k c 
is the crossover scale. In the critical regime, the self-energy correctly approaches the asymptotic 
behavior E(fc) tx k 2 ~ v , and in the short-wavelength regime the behavior is oc fc 2 ' 15-3 ' in D > 3. 
In D — 3, we recover the logarithmic divergence E(fc) oc \n{k/k c ) encountered in perturbation the- 
ory. Our approach yields the crossover scale k c as well as a reasonable estimate for the critical 
exponent n in D — 3. From our scaling function we find for the interaction-induced shift in T c in 
three dimensions, AT C /T C = 1.23an 1//3 , where a is the s-wave scattering length and n is the density, 
in excellent agreement with other approaches. We also discuss the flow of marginal parameters 
in D — 3 and extend our truncation scheme of the renormalization group equations by including 
the six- and eight-point vertex, which yields an improved estimate for the anomalous dimension 
r\ ~ 0.0513. We further calculate the constant lim^o T,(k)/k 2 ~ ri and find good agreement with 
recent Monte-Carlo data. 
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I. INTRODUCTION 

The physics of weakly interacting bosons has seen a 
revival thanks to the improved experimental technique 
of atom traps which allow for a detailed study of Bose- 
Einstein condensation in a controlled environment. Be- 
sides the theoretical effort invested into studies of har- 
monically trapped bosons (for a review, see Ref. 0), this 
also generated renewed interest into the behavior of ho- 
mogeneous Bose gases and led to some new insights 0). 
Many aspects of weakly interacting Bose gases have been 
understood for quite some time now, such as the uni- 
versality class 0(2) of the condensation transition along 
with an accurate knowledge of critical exponents. How- 
ever, the universality of weakly interacting bosons at the 
critical point is not limited to critical exponents since 
a weakly interacting Bose gas behaves universal at all 
length scales larger than the thermal de-Broglie wave- 
length Ath 0- H ■ The origin of this extended universality 
is that at large length scales quantum fluctuations be- 
come unimportant so that this regime of the Bose gas 
is completely described by a classical </> 4 -model p|. In 
this work we present a detailed study of the momentum- 
dependence of the self-energy at the critical point of Bose- 
Einstein condensation at zero frequency, using the func- 
tional renormalization group formalism in the form intro- 
duced by Wetterich Q and by Morris 0. Some results 
of this manuscript were already presented in a brief form 
Here we give a detailed account of the calculation 
and further include an extensive treatment of marginal 
terms. 

In the limit of weak interactions, parametrized by the 



s-wave scattering length a, the self-energy is universal not 
only in the limit of small wave vectors, k — > 0, where it is 
proportional to k 2 ~ ri with a finite anomalous dimension 
■q. It remains universal, in the sense that it can be writ- 
ten in a scaling form independent of a, up to momenta 
which only need to be small compared to AJ^ 1 . At some 
crossover scale k c the zero frequency self-energy leaves 
the anomalous scaling regime and enters the perturbative 
regime, where the self-energy correction to the dispersion 
becomes negligible compared with the bare dispersion 

£ k = Pok 2 , with po = h 2 /2m . (1-1) 

Here m is the bare mass. At finite temperatures and 
D < 4, a perturbative calculation of the self-energy is 
ultraviolet (UV) divergent. This divergence can be ad- 
dressed by using the thermal de-Broglie length as an 
UV cutoff. In D = 3, however, the perturbative regime 
k c < k < 27r/Ath remains non-trivial even in presence of 
a UV cutoff, since additional logarithmic infra-red (IR) 
divergences appear. As discussed in detail by Baym et 
al. y| , the IR divergence cannot be treated in an ad-hoc 
manner by introducing an IR cutoff, since doing so intro- 
duces an artificial additional scale which directly enters 
quantities which should be universal, such as the inter- 
action induced shift of the critical temperature. The IR 
divergence can be removed by a re-summation using a 
variety of standard many-body techniques, e. g. bubblc- 
or ladder-summation and/or self-consistent approaches 
However, these methods are uncontrolled in the crit- 
ical regime where one is faced with a strong-coupling 
problem. Renormalization group (RG) techniques are 
expected to perform better and several authors have ap- 
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plied RG techniques to inve stig ate the IR behavior of 
weakly interacting bosons j^, EH Ell , though no attempt 
was made to calculate the momentum dependence of the 
self-energy. Note that standard field theoretical RG is 
confined to the critical regime k <C k c . In fact, even the 
scale k c cannot be obtained within such an approach. 
To interpolate between the critical and the short wave- 
length regime, functional RG techniques 0,0 are a nat- 
ural choice, since they track the flow of complete vertex 
functions rather than just a small number of coupling 
parameters. Note that the functional RG approach in- 
cludes a priori also terms which are irrelevant according 
to their scaling behavior, but which are important for a 
correct description of the perturbative k c <C fc regime. In 
this work we present in detail a functional RG approach 
to this problem and calculate the universal self-energy to 
leading order in the small parameter an 1 / 3 for all wave 
vectors (n is the boson density). Since the critical behav- 
ior at T c is classical, we will focus only on the self-energy 
E(k, iuj n — 0) at vanishing Matsubara frequency. 

We begin with the general form of the functional RG 
flow equations for the irreducible vertices up to the four- 
point vertex in Sec. [HI where we still retain all frequen- 
cies. Our approach is based on the sharp-cutoff version 
of the functional RG 0HOE1. In SecHTB we turn to 
the effective classical field theory and rewrite the prob- 
lem in a notation appropriate to the classical limit. In 
Sec. lIII Al the relevant and marginal parameters are clas- 
sified in dimensions 3 < D < 4 and their functional flow 
equations are stated in Sec. IIII Bl 

Sec. II VI contains the central part of our work, the cal- 
culation of the self-energy. We first re-examine pertur- 
bation theory and its divergence in D — 3 before we turn 
to the calculation of the self-energy within the functional 
RG formalism. This formalism allows us to calculate the 
scaling function a(x) characterizing the zero frequency 
self-energy at the critical point. We define <r(x) by 

a(x) = (p k 2 c )- 1 mk c x)-i:(0)], (1.2) 

such that the crossover occurs at x = 1. Here S(fc) = 
Ti(k,iuj n — 0) is the exact zero frequency self-energy at 
the critical temperature. Standard field theoretical RG 
can only describe the asymptotic limit x — > 0, where 
a{x) cx x 2 ^ 71 . 

Our approach is based on a truncation of the exact hi- 
erarchy of functional RG flow equations at the four-point 
vertex, i. e. we ignore six-point and higher order vertices. 
While in this approximation the flow of marginal terms is 
not consistently described, the resulting approximation 
for the four-point vertex does in fact include marginal 
and infinitely many irrelevant terms. The resulting flow 
equation for the self-energy can then be solved and we 
derive the complete momentum dependence for the self- 
energy in 3 < D < 4, see Eqs. <|4.51l I4.52|) below, which 
constitute the central result of this work. A numerical 
evaluation of Eqs. (|4.51lH32l in D = 3 is also presented. 
We use a(x) to calculate the shift of the critical tem- 
perature T c of the condensation transition in Sec. In 



D = 3we obtain AT C /T C = 1.23 an 1 / 3 to lowest order in 
an 1 / 3 , in good agreement with recent numerical investi- 
gations 0, 0] and other analytical results 0, 0] (for 
a recent review on this topic see 0)- 

In Sec. lVIl we improve upon the truncation of the flow 
equation and account for the coupling parameters which 
become marginal in D = 3. Three additional parame- 
ters must be taken into account, two of which are associ- 
ated with the linear momentum dependence of the four- 
point vertex and one which describes the momentum- 
independent part of the six-point vertex. We discuss dif- 
ferent truncation schemes of the flow equations and show 
that the inclusion of marginal terms provide an improve- 
ment for the fixed point value of the anomalous dimen- 
sion. 

In Sec. IVIII we calculate <j(x) in D = 3 including the 
marginal terms of the four-point vertex but ignoring irrel- 
evant terms. We demonstrate how this truncation fails in 
the large x regime, where it predicts incorrectly <j(x) cx x 
which in turn would predict AT C /T C cx an 1 / 3 h^an 1 / 3 ). 
On the other hand, in the critical regime irrelevant terms 
only lead to a renormalization of marginal and relevant 
ones and <j{x) is well described by a theory were irrel- 
evant terms are not included. We use this approach to 
express the prefactor A% of the anomalous scaling term, 
cr(x) w A3X 2 ~ V , as a function of r\ for D = 3. Our re- 
sult is in good agreement with recent Monte Carlo results 
4J. Finally, in Scc. lVIlil wc summarize and conclude this 
work. 



II. FUNCTIONAL RG FLOW EQUATIONS FOR 
BOSONS 

Our starting point is a standard effective action de- 
scribing free bosons with a two-particle interaction, 
which is local at the bare level. The action with an UV 
cutoff Ao (to be specified later) is of the form 

S A M,^} = s° Ao {t>,i>} + S^WM , (2-1) 

where the non-interacting part is given by 
S° Ao {$^} = I ©(Ao-|k|) 

JK 

x [-ia; n + e k -/i + £(0,i0)]V5 i W'K . ( 2 - 2 ) 

with <d(x > 0) = 1 and Q(x < 0) = 0. The interaction 
part is given by 

/ e(Ao-|k|)[s Ao (^)-s(o,iO)]^K 

JK 

+ wIkJkJkJ Ki Sk[+k - km 

x Tf o {K[,K' 2 - K 2 , K^k^k^k^k, + . . . , (2.3) 

where the ellipsis denotes three-body and higher order in- 
teractions, which we ignore at the bare level. Here tpx is a 
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complex bosonic field. We use the notation K = (k, iw n ), 

Ik = (PVy 1 Sk,w„' and = /^kk'<L„w„,, where 

(3 is the inverse temperature, V is the volume, and u> n = 
2nnT are bosonic Matsubara frequencies. In Eq. (|2.2|) we 
have included the exact self-energy at vanishing momenta 
and frequencies as a counterterm in the definition of the 
free action. Throughout this work, we shall work at tem- 
peratures T > T c such that the U(l) symmetry is not 
broken. The generating functional of the one-particle ir- 
reducible n-point vertices T A 2n \K' ly . . . , K' n , K n , . . . , K\) 
of the theory with IR cutoff A can be expanded in terms 
of the fields <f>K = H'k) as follows 

1 n f f 

T A {4>,4>} = II / / 5 K' 1 +...+K' n< K„+...+K 1 

y.T { l n \{Kl K s })4> Ki ■ ■ ■ $K>J Kn ■■■<f>K,. (2.4) 

The functional RG flow equations of the first few ir- 
reducible vertices (up to the six-point vertex) for non- 
relativistic fermionic and bosonic many-body systems 
can be found in Ref. ^| • We summarize below the flow 
equations relevant for bosons up to the four-point ver- 
tex. The flow equations for the six-point vertex and some 
terms of the flow of the eight-point vertex are given in 
Appendix [X] 



A. Free energy 



(2) 

where T A (K) is the two-point vertex, defined by 



rf (K) = £ A (AO-£(0,zO) 



(2.6) 



B. Two-point vertex 

The main interest of the work is the flow equation of 
the two point vertex, which we will relate in Sec. lIVI to the 
zero energy scaling function of the self-energy, Eq. (|1.2|1 . 
The flow equation for the two-point vertex is given by 



d A T%\K) 



(2.7) 



A"' 



where 



G A (K) 



<5(A-|k|) 



(2.8) 



iu n -e k + f i--E(0,iO)-T ( l ) {K) _ 
is the cutoff dependent single scale propagator, with sup- 
port only at |k| = A. As is evident from Eq. (|2.7|) . the 
flow of the two point vertex depends on the properties of 
the four-point vertex at finite wave vectors and we thus 
need to derive the momentum-dependent flow of the four- 
point vertex. 



For completeness, we list here the flow equation for the 
free energy T A \ although below we shall not discuss it 
further. The flow equation is given by 



C. Four-point vertex 



x In 



V I 5(A 

IK 



k|) 

iuj n — ek + H — S(0, zO) 
-e k + M -E(0,iO)-rf (K) 



(2.5) 



The flow equation for the irreducible four-point ver- 
tex is responsible for the cross-over from the weak cou- 
pling regime at small scales to the critical regime at large 
scales. The flow is given by 



K 



K 



K 



G A (if)rf(iS-i,^, K- K, K^Kx) 



K 



K'=K 1 +K 2 -K 



G A (K)G A (K')T^ {K'^K'^K 1 , K)vf (K, K' , K 2 , K,) 
[G K (K)G K (K r ) + G A (K)G A (K')]T^(K[,K'; K, K^vf \K' 2 , K; K', K 2 ) 
[G\(K)G\(K') + Ga(K)G a (K')]T^(K 2 ,K'; K, K^vf {K[,K; K', K 2 ) 



K'=K 1 -K' 1 +K 



K>=K r -K'+K 



(2.9) 



where T A is the six-point vertex and 



G A (K) 



9(A < |k| < A ) 
^ n -e k + M -S(0,iO)-rf (K) 



(2.10) 



is the cutoff regularised propagator. The notation 0(A < 
|k| < A ) is shorthand for 6(|k| - A)- 9(|k|-A ). This 
equation is shown graphically in Fig. ^ For the most 
part of this work, we shall in fact ignore the contribution 
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BCS 



of the classical sector by finite Matsubara frequencies al- 
together, since this only leads to corrections in T c of order 
ft 2 . It is convenient to write the single-particle Green's 
function at zero Matsubara frequency in the following 
scaling form 



G A (k,iQ) 



Zi 
PoA 2 



Gi (kA) 



(3.3) 



FIG. 1: Diagrammatic representation of the flow equation for 
the four-point vertex, see Eq. I|2.9[l . We adapt the notation 
BCS, ZS and ZS' for the diagrams from the usual fermionic 
language pj|. even though this notation does not imply a 
physical correspondence. 



from the six-point vertex to this flow and work solely 
with the two- and four-point vertex. To calculate fixed- 
point properties, we will however also include the six- 
point vertex in Sec. IVII 



III. EFFECTIVE CLASSICAL FIELD THEORY 

To discuss the classical critical behavior, it is sufficient 
to retain only the zero Matsubara frequency part of all 
vertices. In principle, all fields with non-zero Matsubara 
frequencies can be integrated out using the flow equa- 
tions given above, which leads to finite renormalizations 
of the parameters appearing in the effective classical the- 
ory ■ The effective UV cutoff of the classical theory 
is determined by the thermal de Broglie wavelength, 



A 2lT 



Ath 



\j2nmT 



2nh 2 
mT 



(3.1) 



In D = 3 the initial value of the four-point vertex can be 
parameterized in terms of the two-body scattering length 



ft, 



r£ (0,0; 0,0) = 



8irh 2 a 



in 



l6irp Q a 



(3.2) 



A. Functional RG flow equations for the rescaled 
vertices 

We assume that the finite renormalizations due to non- 
zero Matsubara frequencies are implicitly taken into ac- 
count via the initial conditions at scale A = Ao of the ef- 
fective classical theory. As discussed in detail in Ref. [j| , 
to calculate the linear shift of the critical temperature 
(see Sec. 01, one can in fact ignore the renormalization 



Gi(q) 



6(K |q| <e l ) 
#/(q) 



(3.4) 



were the minus sign is introduced in Eq. (|3.3|l to arrive 
at the usual definition of the classical Green's function. 
Here / = — ln(A/Ao) is the logarithmic flow parameter, 
and the inverse dimensionless propagator is 



^(q)=Z i q 2 +fp ) (q) + — l - 



[E(0,iO)- M ], (3.5) 



with the dimensionless irreducible two-point vertex 



fP } (q) 



Zl rf (Aq,zO) 



PoA 2 

Zi 
PoA 2 



[E A (Aq,*0)-E(0,tO)]. (3.6) 



The classical wave-function normalization factor Zi is 
given by 



Zr=l- 



af| 2) ( q ) 



dq 2 



(3.7) 



q 2 =0 



Similar to Eq. fity . we define the classical dimension- 
less higher order vertices for n > 2 by 



1 In 7 



■ ,qij 



(k d t) 



n-l 



xA D (»-i)-»»(Z,/po) n ri 2n) ({k 4 =Aq f ,w n< =0}) . 

(3.8) 

For later convenience we have included the numerical fac- 
tor Ko = fl£)/(2Tr) D in the definition of the vertices, 
where flu = 2tt l ) l 2 /T(D /2) is the surface area of the 
D-dimensional unit sphere. The rescaled vertices satisfy 
functional flow equations of the form 



2o. 



2n - D(n - 1) - nr)i - ^ qj ■ V c 



({q,}) 



+rf n) (H}) 



where 



m 



-d,\nZ, 



(3.9) 



(3.10) 



is the flowing anomalous dimension. In particular, the 
rescaled two-point vertex satisfies 



flkf P } (q) = [2 - m - q- VjrHq) + rf'(q) , (3.11) 



(2), 



.(2), 
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where and 

if (q) = / GKqOf^Wjq'.q), (3.12) G ;(q) = ^Lil . (3 .14) 

with The functional flow equation for the dimensionless four- 

/. j£) / point vertex follows from Eq. I|3.9[l . It involves the inho- 

J = J ' ^ 3 ' 13 ^ mogeneity 



J 



f, W (qi,qi;qa,qi) = / G,(q)f f j (q' 1; q 2 , q; q, q 2 , Ql ) - / G,(q)G,(q')f J 4J (qi, q 2 ; q', q)f ; W (q, q'; q 2 , q x ) 



q'=qi-qi+q 



q -qi-q^+q 



q'=qi+q2-q 



[Gl (q)G z (q' ) + G l (q)6, (q')] f | 4) (qi , q'; q, qr )f < 4) (q 2 , q; q' , q 2 ) 
[Gi (q)G z (q' ) + G t (q)6, (q')] f | 4) (q' 2 , q'; q, qj )f[ 4) (qi , q; q' , q 2 ) 

I 



(3.15) 



The flow equations for the six-point and eight-point ver- 
tices also have the form Eq. (|3.9|1 where f and rj 8 - 1 in- 
volve various combinations of the two-, six-, four-, eight - 
and ten-point vertices. All terms entering the inhomo- 
geneity of the six-point vertex and the terms needed 
for a calculation up to second order in the relevant and 
marginal parameters of the inhomogeneity of the eight- 
point can be found in Appendix 1X1 see Eqs. HA2IA3|) . 



with scaling dimension e = 4 — D. The exact flow equa- 
tions of these parameters are 



dm = (2 - 77, )r, + if (0) 



(3.18) 



d m = (4 - D - 2 m ) Ul + r[ 4) (0, 0; 0, 0) . (3.19) 



B. Classification of coupling parameters 



2. Marginal coupling parameter 



Although we are ultimately interested in the flow of 
vertex functions, it is useful to first consider the flow of 
marginal and relevant coupling parameters, since irrele- 
vant parameters become local functions of the relevant 
and marginal ones at the fixed point |20j |. To properly 
organize the flow of irrelevant terms it is thus necessary 
to know the flow of the relevant and marginal parame- 
ters. We first investigate here the relevant and marginal 
terms for 3 < D < 4, ignoring additional marginal terms 
in D = 3. 



1. Relevant coupling parameters 

In 3 < D < 4 there are two relevant coupling parame- 
ters, 



n = lf (0) = lim — i I [S A (Aq,i0)-E(0,i0)] , (3.16) 

q^O p A z 

with scaling dimension +2, and 

ui = f[ 4) (0,0;0,0) , (3.17) 



For 3 < D < 4 the only marginal parameter is the 
wave- function renormalization Z\. The exact flow equa- 
tion is 



diZi = -riiZi , 
where the flowing anomalous dimension is 



m 



9ff (q) 



and 



<->a - 



d D q 



*(|q|-i) 



(3.20) 



(3.21) 



(3.22) 



denotes the integral over the unit sphere. The surface 
area H,d of the unit sphere was defined below Eq. (|3.8|l . 
In D = 3, additional marginal coupling parameters ap- 
pear which we ignore in this section. We will discuss 
them in Sec. IVII 
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3. The Wilson-Fisher fixed point close to D = 4 

It is instructive to take a closer look at the perturba- 
tive one-loop RG flow equations for the relevant coupling 
parameters. If we retain from the two-point vertex only 
the relevant and marginal part, we may approximate 



so that 



Gi(q) 



i?i(q)~n+q 2 , 
9(1 < |q| < e l ) 



G;(q) 



n + q^ 

g(|q|-i) 

n + q 2 



(3.23) 



(3.24) 



(3.25) 



If in addition we retain only the relevant part of the 
four-point vertex, we may replace it by a momentum- 
independent constant, 



f [ 4) (qi , q' 2 ; q 2 , qi ) « f \*> (0, 0; 0, 0) = m 



( 4 )/ 



(3.26) 



Within this approximation we obtain from Eq. I|3.12|) . 



f(2) 
1 I 



(0) 



Ul 



1 + n 



and Eq. (|3.15|l reduces to 
f| 4) (0,0;0,0) 



2(l + r ; ) 2 



(3.27) 



(3.28) 



Because the anomalous dimension is related to the mo- 
mentum dependence of the four-point vertex which is 
irrelevant in D > 3, we may set rji ss within this 
approximation. Furthermore, on the critical trajectory 
ri = 0(ui) can also be ignored in Eq. 13.28fl as long as 
ui remains small compared with unity. Then we obtain, 
with e = 4 - D, 



Q ° 2 

diui = eut - -u t , 



din = 2r/ + ui 



(3.29) 
(3.30) 



Eq. I|3.29|l is easily solved. The solution can be written 
in the form 



Hi 



1 



1 



1 



where 



e{l-l c ) _|_ l e e(Ic-0 + 1 



2 

u* = -e 
5 



(3.31) 



(3.32) 



is the value of ui at the RG fixed point, and the logarith- 
mic crossover scale is 



Z c = iln(^-lUiln(^ 
e \u I e \uo 



(3.33) 



U* 







FIG. 2: Typical flow of the coupling ui as a function of I for 
wo <K w„ and D = 3. 



Throughout this work we assume that uq <C it*, corre- 
sponding to weak bare interactions between the bosons. 
Note that the right-hand side of Eq. Ij3.31|l is expressed 
in terms of the Fermi function. If we think of e as the 
inverse temperature and l c as the chemical potential, the 
qualitative behavior of ui/u* is clear: within a narrow 
(on the scale l c ) interval of width 1/e centered at I = l c 
the ratio ui/u* raises from the small value e _e ^ c_ '' to a 
value close to unity, see Fig. [21 Given the solution ui in 
Eq. H3.31fl . the RG equation (|3.30|l for r; is easily solved, 



n 



,2/ 



TO 



e 2e 

5 ~ Wo 



dl'- 



-21' 



,c(l'-l c ) _)_ I 



(3.34) 



To obtain a fixed point at I — > oo, the intial value r$ has 
to be fine tuned such that 



2c 



dl'- 



-21' 



4i'-ic) + i 



Then we may write 



n = + u» 



dl'- 



-2(1' -l) 



where 



Y 



■ e(l'-lc) -)- 1 



(3.35) 



(3.36) 



(3.37) 



In Sec. IIVI we show that, up to a numerical factor of the 
order of unity, the momentum scale k c = Aoe _ic associ- 
ated with the logarithmic scale factor l c , i. e. 



k c = A 





-1/6 




It* 


,A (^) 


_M 







l/e 



(3.38) 



can be identified with the crossover scale where the criti- 
cal /c 2_I, -form of the energy dispersion begins to emerge. 
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IV. THE SELF-ENERGY AT THE CRITICAL 
POINT OF BOSE-EINSTEIN CONDENSATION 




A. Second order perturbation theory 



To begin with, let us attempt to calculate the self- 
energy by means of straightforward second order pertur- 
bation theory [|| , which yields in a continuum model with 
UV cutoff A 



E(k) - E(0) = -t\t 2 J -^JXO(P) 



(2tt)' 

e(A -|p + k|) e(A 



(4.1) 



e p+k 



where T = Tj^(0, 0; 0, 0) is the bare vertex, and 



Xo(p)=T 
With e 



d D P ' e(Ao-|p'|)e(A -|p' + p|) 



PoP we obtain for Ao — * oo, 
T 

Xo(p) = K D K D 

PoP 



(4.2) 
(4.3) 



where 



K' = 



dx[x(l-x)] {D - i)/2 / dy 



./(, (y 2 



,D-1 



ya-o v^r (f - 1) fCD-2) 



r(2=l) sin(§(Z>-2) 



(4.4) 



Note that K' D is finite for 2 < D < 4; in particular, K' 3 = 
7r 2 /4. Substituting Eq. I|4.3j) into Eq. (|4.2j) and taking 
the limit Aq oo, we can scale out the /c-dependcncc 
and obtain for D < 3 < 4, using the definitions uo = 
ifuTAg e po 2 T and fc c = Ao(u /w») 1/e , 



PoMW-ECO)] =B D kfk 2 



(4.5) 



where after an integration by parts the coefficient 
can be cast into the following form, 



B D = 



3ui 



A{D - 3) 



K'l 



D 1 D 



(4.6) 



D-l 



n, 



ddismd) - 2 I dxx 

n Jo 

1 + 2x cos $ 



2D-6 



da; V 1 + 2x cos i9 + x 2 



(4.7) 



In Fig. |3| we show a numerical evaluation of Bp for 3 < 
D < 4. With 7 3 = 1 and K' a = vr 2 /4 we obtain for D -> 3 



2, ,2 



D-3 



3tHu 
16 



(4.8) 



0.5 



3.4 3.6 



3.8 D 4 



FIG. 3: Graph of the coefficient Bd defined in Eqs. I|4.6^ and 
114.711 in dimensions 3 < D < 4. 



Perturbation theory is valid as long as the correction l|4.5|) 

is small compared with pok 2 , i. e. for k > k c B^ 2e . Due to 
the IR divergence of Br> for D — > 3, close to three dimen- 
sions perturbation theory breaks down even outside the 
critical regime. In a IR regularized approach with a cut- 
off £, the divergence in Eq. (|4.5I) is replaced by a ln(fc£) 
behavior |5j. Below we shall see that the functional RG 
approach shows a ln(fc/fc c ) behavior, as expected for a 
theory where the divergence due to density fluctuations 
is screened within a non-perturbative treatment. 



B. Functional RG calculation of the self-energy 

We now derive the behavior of the self-energy micro- 
scopically using the functional RG equations for the ef- 
fective classical field theory. 



1, The 



u\ -truncation for the inhomogeneity 



To begin with, we need the RG flow of the four- 
point vertex. In the simplest approximation, we expand 

Fj ; (oiiQa! 12, qi) in powers of ui, which should be accu- 
rate as long as the renormalized coupling remains small. 
To leading order we simply ignore the six-point vertex 
and replace the four-point vertices on the right-hand side 
of Eq. (|3.15|) by their relevant part ui. In this approxi- 
mation, 



r| 4) (qi,q 2 ;q2,qi) ~ -uf ^x/(|qi + q2 



i 

i v (.qi,q 2 ;q2,qij » ~ u t[2 : 

*i(ln-qil) + *j(|qi-qal)l , (4.9) 



where the generalized susceptibility is given by 



8 



XM) = 2 / G(|q'|)G,(|q / + q|) 

Jq' 



2 /9(1< |q' + q| < e l ) 



l + r/\ |q' + q| 2 +r ; / q' 

2 ^_D-l /"^ In/ _. _ _ a \D-2 



1 + ri Q D j 
6(0 < q 2 + 2gcos?J < e 21 - 1) 
1 + n + <Z 2 + 2(7COSi? 



(4.10) 



where in the second and third line we have used the ap- 
proximations (|3.24|1 and l|3.25fl . taking into account only 
the relevant and marginal part of the two-point vertex. 



2. Generalized susceptibility 

It is instructive to examine the behavior of the gener- 
alized susceptibility xi (q) ■ For e) — 1 3> 2 the asymptotic 
behavior for small and large q can be easily obtained an- 
alytically, 



The important point is that the leading correction for 
small q is linear in q: expanding Eq. I|4.10|l in powers of 

q, 



Xi(q)=Xi(0)+Xi(0)q + O(q 2 ) , (4.12) 



we find 



X'i(0) = ~ 



Xi(0) = 



(l + n) 2 



{l+n) 2 ' 



4S 1 ! 



(l + n) 



-S 



where 



and 



So = (S(q' ■ q)> Q/ = 



So 



^ = (e(q'.q)q'.q)., £) _ | 

In D = 3 we obtain with So — 1/2 and Si = 1/4 
1 



(4.13) 



(4.14) 



(4.15) 



(4.16) 



x'M = - 



(1 + r;) 2 \l + r, 



' -J] , D = 3. (4.17) 



For later reference, let us give the exact function Xl(q) 
in D = 3, where the angular integration in Eq. (|4.10|) 



can easily be performed analytically. The result can be 
written as 



. , , Q(g2 - gi) . 



2(1 + rj)? 



1 + n + q 2 + 2qq 2 



l + n+q 2 + 2qqi 



where 



q\ 



92 = 



-1 if g > 2 

-q/2 if q < 2 ' 

1 if e l - 1 > q 



S^-l if' 

2q 2 



Kg 



(4.18) 

(4.19) 
(4.20) 



In particular, for q < min{2, e l — 1} we have 



xi(q) 



1 



2{l + n) q 



In 



{i + q) 2 + n 



l + n 



(4.21) 



Thus, for small momenta our approximation l|4.9|l yields 
in D = 3, 

r} 4) (q'i,q^;q2,qi) ' 



2(l + r,) 2 



1 



(l + n) 2 \l + n 2 
|qi-qil + |qi-q' 2 | 



2 |qi + q2| 



+ 0(q 2 ) 



(4.22) 



The term linear in q generates marginal parameters, as 
mentioned in Sec lIII B~2l when we iterate the RG. The 
linear term exists for all Z?, but the corresponding cou- 
pling parameters are irrelevant in D > 3. Because these 
coupling parameters are not consistently taken into ac- 
count in the u 2 -truncation given in Eq. 14.91) . we cannot 
expect that this truncation gives numerically accurate 
results in the critical regime close to three dimensions. 
On the other hand, in the short-wavelength regime (and 
also for D — 3 ^> 77), our u 2 -truncation 14. 9f) is suffi- 
cient. In this section we shall therefore proceed with this 
approximation, which produces well-defined results even 
in D = 3 provided the anomalous exponent 77 is calcu- 
lated self-consistently by solving an integral equation, see 
Eq. (|4.29() below. In Sec. I VII we shall improve on this ap- 
proximation by explicitly including all marginal coupling 
parameters in the critical regime in D — 3. 



Explicit expression for the four-point vertex 



With Eq. (|4.9() as an approximation for r\ , we are 
now in a position to calculate the flow of the four-point 
vertex. Let us rewrite Eq. I|4.9fl in the following way, 



r| 4) (q , i,q 2 ;q2,qi) 



-5u 2 
2(1 + nf 



r, v (qi,q 2 ;q2,qi) > 



(4.23) 
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where 



r ; (qi,q 2 ;q2,qi) = -u x 



with 



(mi). 



+x| mi) 



qi-q'il)+x| mi) (|qi-q2l)] , (4-24) 



describing the momentum-dependent part. The super- 
script (mi) indicates that these terms are marginal or 
qi + q2 1 ) irrelevant. In this approximation the four-point vertex 

takes the form 

f ; (4) (q , 1 ,q^;q 2 ,qi) «^+f ; (4mi) (q 1 ,qi;q2,qi) , (4.26) 



with the irrelevant (and in D — 3 also marginal) parts 
(4.25) given by 



fS(4mi)/ / / \ 

ij (qi,q 2 ;q2,qi) 



^ e e(I-n-Jj;^fW( e -Wql )e -Wqi;e-Wq !ie -( , - , \) 



dl > e e{l-l')-2Jl,dT^ u 2 



k^V^lq! + q 2 |) + X^ie-^hi ~ qil) + X^V^'lqi - <&\) 



.(mi)/ — 1\ \\ , -(mi)/ —t\ f\\ , .(mi)/ —t\ / 

Xi- t (e |qi + q2|)+4_ t ( e Iqi -qil) + X/_/(e |qi-q 2 



.(mi) 



(mi) , 



(4.27) 



C. Anomalous dimension 

To obtain the anomalous dimension r\i = 
dtf\ci)/dq 2 \ q 2 =a we calculate f f } (q) via E q. fTTfy 
With the approximation Eqs. (|4.26l I4.27|) for the 
four-point vertex we find 



fp } (q) = 



1 



/?/(:| ^f^(q,q';q',q; 



q' 



ui 



T £ 



\ + n 2{l + n )J 

x (^(e-'dq'+ql))-, • (4.28) 

From Eq. 1|3.21|) we then obtain an integral equation for 
the flowing anomalous dimension 



m 



with 
K(l,t) 



dtK(l,t)u 2 _ t e- 2 tt-t dTr >T , (4.29) 

^^^(^(e-^lq' + ql))).,!^ 
p-l)e- (D - 3)t xU(e- f ) 

(4.30) 



4L»(l + r/ 

+e-^x'Ue- t ) 



where xi( x ) = dxi(x)/dx and x'ii x ) — d 2 \i{x) / dx 2 . 
Note that on the right-hand side we have replaced 

Xi-til) ~ * Xl-t(q), because the constant part xi-t(0) 
does not contribute to the derivative. The integral equa- 
tion (|4.29() together with the one-loop flow equations for 
the relevant coupling parameters, 

d in - (2 - m ) n + -HL_ (4.31) 

1 + n 




FIG. 4: Self-consistent RG flows of ui, r; and r/i at the critical 
point in D — 3 (dashed lines), obtained numerically from the 
coupled intego-differential Eqs. I|4. 2914. 3114. 32^ . For compari- 
son, the Fermi-function approximation Eq. I|4.33|l for the flows 
is also shown (full lines) , with fixed point values matching the 
numerical solution. 



dim = (e - 2i]i)ui - - 



2(i + n) 2 ' 



(4.32) 



form a system of three coupled integro-differential equa- 
tions for the three unknown functions ri, uj, and r/i. A 
numerical solution of these equations for D = 3 is shown 
in Fig.H 

To leading order in e (i. e. close to four dimensions) 
Eqs. (|4. 2914. 3114. 32fl decouple and one recovers the fixed 
point values it* = 2e/5, = — e/5 from Eq. I|3.37|l . 
To this order it is consistent to ignore rji and r; on 
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the right-hand side of Eq. i|4.32[l . so that it reduces to 
Eq. i|3.29[) . with the Fermi function solution for it; given 
in Eq. H3.31|l . From the numerical solution of Eqs. 1)4.291 - 
14.321) for D = 3 we find that even if r\ and rji are taken 
into account, to a good approximation the qualitative be- 
havior of it; is still well described by a Fermi function if 
we adjust the fixed point value it* to the one obtained 
from the numerical solution. Furthermore, Fig. 0] shows 
that also r\\ and r; roughly follow the same functional 
form and we therefore approximate 



m 

V 



11 

r* 



1 



(4.33) 



with the fixed point values taken from the numerical so- 
lution and where r\ = lim;^,^ 77; is the anomalous dimen- 
sion at the fixed point. Although the functional form 
(|4.33l) for rji is qualitatively incorrect for I <C l c (where 
i]i oc uf, whereas Eq. (|4.33|) predicts r\i cx u{) we may use 
Eq. H4.33f) as a zeroth approximation on the right-hand 
side of the integral equation (|4.29() . The point is that for 
I <C l c the right-hand side of Eq. 14.29J1 is not sensitive 
to the precise value of r\i in this regime, and the first it- 
eration of the integral equation yields an accurate result 
even for small I. Let us now calculate the fixed point 
values. We obtain from Eq. 



V = K I dtK (oo, t)e 
lo 



-2r)t 



(4.34) 



where we use the notation r\ = 77*. The function K (00, t) 
is obtained from K(l, t) defined in Eq. (|4.30|) by replacing 
the functions Xi-t(q) on the right-hand side by 



X*(q) = hm xi(q) 



6(|q' 



1) 



In particular, in D 



3 we have 



(4.35) 



x*(q) 



1 



6(2 - q) In 



2(l + r*)g 

+G(q - 2) In 



(l + q) 2 +r, 



1 + r* 



(1 + 9) 2 



(1 -q) 2 + 7% 



(4.36) 



Note that for large t the kernel K (00, t) vanishes as 
e _ ( D_3 ^*, so that in D = 3 it is crucial to retain the 
77-dependence on the right-hand side of Eq. I|4.34|) . This 
is closely related to the appearance of marginal terms of 
the four-point vertex in D = 3. The self-consistent fixed 
point values can be obtained by using the fact that, at 
the fixed point, we can express and r* as a function 
of T), 



2 1-2^ 
_ 3 4 - 377 ' 
10 (1-2t ? )(2-t ? ) 2 
9 (4 - 3t]) 2 



(4.37a) 
(4.37b) 



7% p 


a -0.143, 




d 0.232 , 


T] P 


s 0.104. 



Inserting these expressions for it* and r* into Eq. I|4.34(l . 
we obtain a self-consistent equation for r\ which can be 
solved numerically. The numerical solution gives 

(4.38a) 
(4.38b) 
(4.38c) 

The above value for 77 is approximately three times as 
large as the generally accepted valu e n « 0.038 obtained 
by several different techniques [22,|2j|]. However, given 
the simplicity of our truncation, it is quite satisfactory 
that our estimate H4.38c|l for r\ has the correct order of 
magnitude. In any case, our simple u 2 -truncation cer- 
tainly gives a much better estimate for 77 than the self- 
consistent two-loop calculation employed by Baym et al. 
|5j, which predicts a value of 0.5 for 77. In Sec I VII we shall 
further improve on the it 2 -truncation employed here by 
explicitly taking the RG flow of marginal coupling pa- 
rameters into account. 

For completeness, we show how to obtain the standard 
result for 77 to second order in an e-cxpansion. We can 
rewrite Eq. I|4.34|l as 



Sui 



4£>(l+r*) J 



d\\- e+271 



(D-1) X {(X)+X X ':(X) 



4D(l + r») 



2(1-77) / dAA- £+2 V*(A)+x;(l) • 


(4.39) 



For small e = 4 — D we obtain to leading order 
^ = -^[2(*(l)-*(0)) + *'(!)] » 



(4.40) 



where the function x(q) is obtained from X*{q) by simply 
setting r* — > 0. Evaluating the integrals in D = 4 (see 
also Ref. [13), 



X(0) = 1 , 

4 _ V3 
3 ~ V 



X(l) 



■>m 4 2^3 

S 7T 



(4.41a) 
(4.41b) 

(4.41c) 



and using it* = 2e/5, we finally obtain 77 = e 2 /50 + O(e 3 ), 
in agreement with the field theoretical result [2l| . 



D. The scaling function 

Let us now calculate the dimensionless scaling function 
<Ji{x) which we define by 



07 (x) 



-w-Uz^rPieV-Ux) 



(4.42) 



Since the two-point vertex (q) and the inhomogeneity 
T\ ' (q) depend only on q — |q| , we shall in this subsection 
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use scalar arguments for these functions to simplify the 
notation. For I — ► oo the definition i|4.42H for <ti(x) yields 
the universal scaling function defined in Eq. 1)1.2(1 . 

lim ai(x) = a(x) = (pohgy^kex) - E(0)] . (4.43) 

I — >oo 

We first perform some exact manipulations. The func- 
tional flow equation for f[ 2 \q) given by Eq. I|3.11fl can 
be transformed into an integral equation 

+ f dte^-ti-^ff^e-tq) 
Jo 



a 2i-/' drr, T 



. (4.44) 



To describe a critical system, we choose the initial value 
M=o(<7) = r o to be momentum independent such that for 
I — > oo the relevant coupling 



+ / dz'e- 2i ' + /o'<^r, ( , 2) (o) 



(4.45) 

has a finite limit r* = lim^oo r;, just has we have done 
in Eq. (|3.34|) . This is guaranteed if the initial value ro is 
chosen such that l24l 



ro 

Defining 



r ( 2mi ) ((z) = rW(9) _ rr(0)j 



(2), 



,(2), 



(4.46) 



(4.47) 



where the superscript mi indicates that these function 
contain only marginal and irrelevant parameters, we ob- 
tain 



d l> e -W'-lc)+fl'dT VTt (2^( e l'-h x ) (448) 



The first term on the right hand side vanishes for I — > oo 
because by construction lim;_ ! . 00 ri = r* is finite on the 
critical surface. Taking the limit I — > oo we thus obtain 
(after renaming I' — > /) 

/>oo 

*(x)= dle-^ l - l ^o dr ^tf m ^{e l -^x) . (4.49) 



1. Truncation of the flow equations 

So far, no approximation has been made. We now 
approximate the function r| 2mi ) (g) n the right-hand side 



of Eq. H4.49fl by the leading term in the expansion in 
powers of ui, see Eq. (|4.28[) . 

-A(2mi), s 3 f ei-2 f, ! , d,Tri T 2 

xixi-tie-'W + ^-Xi-tie-'))^, 
2(1 + n) J 1 

x(x,(e- ( '- n |q' + q|) -Xp(e- ( '- n ))-, . (4.50) 

All quantities on the right-hand side of Eq. (|4.49|) are now 
known, so that we may calculate the crossover function 
by performing the four-dimensional integration (two an- 
gular integrations and two integrations over the scale pa- 
rameters I and I') numerically. There are no divergences 
even in D = 3 provided the integral equation l|4.29|l for 
r/i is solved self-consistently. 



2. Results for the scaling function 

To obtain an analytic approximation for the scaling 
function a(x), we adopt again the successful strategy 
used in the solution of the integral equation (|4.29(l for 
n: we substitute the Fermi function ansatz (|4.33l) for the 
flowing anomalous dimension n T on the right-hand side 
of Eqs. (|4.49|) and Q4.5Up . After some transformations of 
the integration variables we then obtain 



d y-< 



y 



-3+2c 



+ yeU-2 V /e 



F(x,y;r),l c ) 



(4.51) 



dz- 



a(x) 
with 

F(x,y;r),l c ) — . . 

la [x- + (y/zY] n/£ [l + r ic+ln ^_ 

x(xi c+ i n |(^) -Xl e +l**(\z$ + V$'\)) v > ( 4 -52) 

where q" is an arbitrary unit vector. Note that by as- 
sumption x -C e lc so that we may replace xe~ lc — ► in 
the lower limit of the y-integral in Eq. (|4.51|) . Consider 
first the regime x <C 1. Using the fact that for I ^S> l c we 
may approximat e n — > and ±M) — * lini/^oo ±i (?) = 
X*{q) ( see Eq. (|4.35|) ~). we may replace in this regime 
Xi c +lnf (?) -> X*(q) and r ic+ini -> r * in Eq. I4.52|l. 
Then we obtain 



3?; 2 r c 



y 



-3+2e 



:F*(x,y;ri) 



2 J ~ a [x £ + ^p- 2 ^ 

Using D-dimensional spherical coordinates, the function 
F*(x, y; rj) = lim/ c _ +00 F{x, y; n, l c ) can be written as 



F*(x,y;v) = 



d Z r Mf-wr-' 



{l + r*)fl D Jo Jo [x e + (y/z) 



X X*(z) - X*Wz 2 + 2zycos$ + y 2 ) . (4.54) 
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From Eqs. (|4.53(l and (|4.54(l it is now straightforward to 
obtain the asymptotic behavior of <j(x) for small x, 



a(x) 



Adx 2 ^ 



(4.55) 



2(1 + r.) Q D J J Jo 



x(sin^) D " 2 x*(z) - X*W* 2 + 2zy cosi? + y 2 ) .(4.56) 



In D = 3 we find numerically 



A, 



1.17. 



(4.57) 



On the other hand, for large x the flow parameter in 
Xlc+hx ^ (<?) an d r ; c +in is typically small compared with 
Z c . To take the effect of the flowing r; approximately into 
account, we replace in this regime r/ c+ i n _«_ — > ro and 
Xz c +lnf (?) -» Xo(q), where xo(g ) is ob tained from x*(o) 
by replacing r* — > ro . Using Eq. (|3.35|) one can show that 
on the critical surface and for small Uq, the initial value of 
ro is approximately given by ro « — 1*0/ (2 — e). With this 
approximation we obtain in the regime [2(D — 3)] _1 <C 
In a;, 



I? 



3u 2 



2(1+ ro) «Z5 



- 3+2e ' dzz 
'0 



l-e 



(4.58) 



dO 



x(sini?) I) - 2 X o(z) ~ Xo( V^ 2 + 2zy cos 1? + y 2 ) , (4.59) 

and for 1 -c lnx <C [2(D — 3)] _1 (which includes the limit 
D->3), 



a{x) 



where 



B' 



3u 2 



2(1 + r ) Jo 



dzXo(z) 



^3 , 



3^ 2 



24 



(4.60) 



(4.61) 



and on the right-hand side we have set ro 



merically we find B$ 
self-energy 



0. Nu- 

0.0319. Note that the physical 



E(fc) - E(0) = p Q k 2 c a(k/k c ) 



(4.62) 



is independent of u 2 in the regime k 3> k c : keeping in 
mind that k c = Ao(uo/w*) 1 / e , it is obvious that the factor 
of it 2 in Eqs. (|4.59|) and 1)4. 61[) combines with a suitable 
power of k c so that the self-energy is proportional to U 2 ,. 
A numerical evaluation of a(x) is shown if Fig. [3] 

To summarize, using the functional RG, we have 
derived in this section the momentum-dependent self- 
energy of weakly interacting bosons, covering the en- 
tire range from the critical regime up to momenta of 




FIG. 5: Numerical evaluation of the scaling function a(x) of 
the self-energy in the limit a — * and D = 3. Dashed lines 
are the asymptotic forms predicted by Eqs. 14.5514.61)1 , 



the order of the inverse thermal de Broglie wavelength. 
While in the critical regime k <C k c it is simpler to ob- 
tain the asymptotic long- wavelength behavior of the self- 
energy within the field theoretical RG, the functional RG 
approach adopted in this work can describe the entire 
crossover from the critical regime k <C k c to the short- 
wavelength regime k c <C k < 2-w/Xth- Note that the be- 
havior of the self-energy in the short-wavelength regime 
is determined by irrelevant terms which are simply dis- 
carded in the field theoretical RG. 



V. INTERACTION INDUCED T C -SHIFT 

Baym et al. Q have shown that to lowest order in 
the scattering length a, the interaction induced shift for 
the critical temperature T c in D — 3 can be calculated 
from the k-dependence of the zero frequency self-energy 
E(k, 0). In a straightforward generalization of their result 
to arbitrary dimensions, we can write the contribution 
from classical modes to the shift in T c as 



A /k c 



211 



D 



D-2 



dxx 



nD((D/2) \A Q/ 

d-3 (^-^(M-sCo)] 



.T2 + ( Po fc2)-i[E(fc c x)-E(0)] 



(5.1) 



Using k c /A Q = e~ Lc = (up/it*) 1 / 6 , and substituting our 
scaling function <r(x) as defined in Eq. (|1.2(l . this can also 
be written as 



T c 



2Qp ( uo 
ttD((D/2) \u* 



D-2 
4-D 



A /fc c 



dxx 



D-3 



x 2 + <j{x) 



(5.2) 
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The shift in T c is dominated by classical fluctuations 
as long as the UV cutoff A /fc c can be removed in 
Eq. (|5.2(1 . Keeping in mind that according to Eq. (|4.58l) 
a(x) ~ Bdx 2 ( d ~ 3 ^ for large x, we see that the integral 
(|5.1|) is UV convergent as long as < 2, i. c. D < 10/3. 
For D > 10/3 the value of the integral in Eq. I|5.1|l de- 
pends on the UV cutoff, such that the lowest order con- 
tribution to the shift in T c is proportional to itg, with 
logarithmic corrections (oc UqItiuq) in D = 10/3. Hence 
for D > 10/3 the shift in the critical temperature cannot 
be obtained from an effective classical field theory, while 
for D < 10/3, the cutoff only leads to higher order cor- 
rections (in D = 3 one finds Uq and Uq\tl(uq) corrections 
(3). Writing 



T c 

we obtain for D < 10/3 
2Q, D 



= JdUq 



(5.3) 



Jd = 



ttD((D/2) 



dxx 



D-3 



a(x) 



a{x) 



(5.4) 

Note that Jd implicitly depends on uq via a{x). How- 
ever, for uq — > the coefficient Jd approaches a finite 
limit independent of u as long as D < 10/3. 
In D = 3 one usually writes 



c\an 



1/3 



(5.5) 



Keeping in mind that 

mo = l6ir- 1 [({3/2)]- 1/3 an 1/3 for D = 3 , (5.6) 

we have in this case c x = 167r _1 [C(3/2)] _1 / 3 J 3 . Within 
our approximation we find c\ — 1.23, in rather good 
agreement with the most accurate numerical investiga- 
tions which give ci = 1.30±0.02 14] and c x = 1.29±0.05 
|15| . Variational perturbation techniques give similar re- 
sults, ci = 1.23 ±0.12 16], and c x = 1.27 ± 0.11 PH. 
More results from other approaches are summarized in a 
recent review |2J . There is also experimental data from a 
4 He-V yco r system confirming the linear scaling of T c with 
an 1 / 3 |25| . The reported value c± « 4.66 seems however 
rather large which may be attributed to uncertainties re- 
garding the precise value of a in the experiment, see the 
discussion in Ref. 1151. 



VI. IMPROVED DESCRIPTION OF THE 
CRITICAL REGIME IN D = 3: INCLUDING 
MARGINAL TERMS 

Within the w 2 -truncation employed so far, we have ne- 
glected the six-point vertex and approximated the inho- 
mogeneity in the flow equation of the four-point vertex 
by 



r| J (qi,q' 2 ;q2,qi) w -Ui 

+Mhi - q'il) + 



2»(|qi 



q 2 |) 

XidQi-Qal) » 



(6.1) 



see Eq. (|4.9|) . This approximation amounts to a trunca- 
tion of the exact hierarchy of flow equations where on the 
right-hand side of Eq. H3.15(l we simply replace 



f ; (4) (qi,q2;q2,qi) 



(6.2) 



A formal justification for this procedure can only be given 
for D > 3: for weak interactions, irrelevant coupling 
parameters can be expanded in powers of relevant and 
marginal ones fl2l l2fj| . For 3 < D < 4 the only part of 
the four-point vertex which is not irrelevant is its con- 
stant part it;. However, in D = 3 there are two ad- 
ditional marginal parameters related to the momentum 
dependence of the four-point vertex. Hence, in this case 
Eq. H6.1fl is not consistent. For small momenta the flow- 
ing four-point vertex has the form 



f; (4) (qi,q' 2 ;q2,qi) = «z 



-a-i (|qi 

-1/ (qi,q 2 ;q2,qi 



q'J + lqi -q^|) + 6 ; |qi+q 2 | 

(6.3) 



where the irrelevant part f ^ (q' 1; q^; q 2 , qi) vanishes at 
least quadratically if all momenta become small. Below, 
we will denote the diagrams entering the renormaliza- 
tion of the four-point vertex using the usual fermionic 
language (i. e. BCS, ZS and ZS', see Fig.^l. We empha- 
size that this terminology simply describes the topology 
of the diagrams, not the physical phenomena arising from 
the diagrams. Note that the four-point vertex must be 
symmetric under the exchange qi <-> q 2 and q^ <-> q^, 
so that both "zero sound channels" are characterized by 
the same parameter ai. The parameter bi characterizes 
the BCS channel. 

The exact flow equations for the coupling parameters 
ai and bi are 



di ai 
dih 



(3 

(3- 



D 
D 



2tji)aj + rf 



1 1 i 



(6.4) 
(6.5) 



where the coefficients f f and are defined via the ex- 
pansion of the function r[ 4 ^ (q^ , q^ ; q 2 , qi ) in Eq. (|3.15|) 
for small momenta, 



r1 4) (qi,q 2 ;q 2 ,qi) = Tj 4J (0, 0; 0, 0) 

+ ffdcu-qil + lqi-qil) 

q 2 |+0(q2). (6.6) 



Itlqi 



Even if the initial vertex fj=Q (qi , q 2 ; q 2 , qi ) is 
momentum-independent (corresponding to the initial 
conditions do = &o = 0) , finite values of these parameters 
are generated as we iterate the RG. Because in D = 3 
these are marginal, they cannot be ignored and are pos- 
sibly the source of logarithmic corrections. 

Besides ai and bi, there is a third marginal parameter 
in D — 3, the momentum-independent part of the six- 
point vertex, 



vi =f ; (6) (0,0,0; 0,0,0) , 



(6.7) 
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see Eq. (|3.9|) . It satisfies the exact flow equation 

8 lVl = (6 -2D- 3 m ) Vl + f | 6) (0, 0, 0; 0, 0, 0) . (6.8) 

Again, the RG flow generates a finite vi even if initially 
wo = 0. At the Wilson-Fisher fixed point in D = 3 the 
renormalized u* = lirn/_>oo v i is of order unity. 

To arrive at the RG equations, we adopt the following 
approximation scheme: We expand r[ n ^ to some fixed 
power in the marginal parameters ai,bi,vi, and the rel- 
evant parameter ui while keeping higher orders in the 
relevant parameter r;. The truncation is based on the 
hope that, although ai, 6j, vi and ui are not paramet- 
rically small, these parameters still remain numerically 
small. The fixed point values for these parameters which 
we obtain within this scheme in Sec. IVID II are indeed 
much smaller than unity, which gives some a posteriori 
justification to our approach. To implement the scheme, 
it is useful to transform Eq. I|3.9[) into an integral equa- 
tion, 

fj rl) ({qj) = e^-^+^-^^f^ae-'q,}) 

+ f dte [(2 - D)n+D]t - n ^ dT ^t^ t {{e-\i\) . (6.9) 
Jo 

We use these retarded expressions in the calculations of 
the inhomogeneities T) m \ Since the right hand side in 
Eq. H6.9fl depends on vertices which also have the form 
1)6.9(1 . one quickly arrives at rather complex expressions 
which we need to truncate. We will only include terms 
up to second order in ai,bi,vi, and u/. However, even 
with this approximation the calculation becomes rather 
involved. To keep the calculations tractable, we truncate 
further: To evaluate the flow of parameters associated 
with the n-point vertex, we ignore irrelevant contribu- 
tions arising from vertices of order m with m > n + 2. 
Thus, when we calculate the flow of parameters associ- 
ated with the n-point vertex and encounter a vertex f 



in the calculation with m > n + 2, we only keep its rel- 
evant and marginal part instead of employing Eq. I(6.9|l . 
Effectively, we neglect in this way infinitely many irrel- 
evant parameters related to higher order vertices. The 
expectation is, that these irrelevant contributions are 
numerically small compared to the contribtutions orig- 
inating from irrelevant terms of lower order vertices. In 
concrete terms, our truncation is as follows: When we 
calculate the flow of vi , we will consider irrelevant contri- 
butions from f but ignore contributions from fj ^ (in 
this case, there are in fact no second order terms arising 
from inhomogeneities T) with n > 10 so that we really 
have all second order terms included in our flow of vi). 
Similarly, to calculate the flow of the four-point vertex 
parameters a/ , bi , and Ui , we ignore contributions arising 
from f; . This leads to the absence of vf terms in the 
RG equations for a/ , 6; , and it; , the T) terms only enter 
indirectly via the flow of Vi . To calculate the flow of r; we 
ignore irrelevant contributions from f > for n > 6 which 
again only enter via their contributions to the flow of a; , 
bi, and ui. 



A. Truncating the flow equation for the six-point 
vertex 

We begin with the renormalization of vi, the momen- 
tum independent part of f j , using the flow equations 
for the six-point vertex from Appendix^ Eq. (|A2|) . Note 
that, working consistently to second order, the last two 
terms in Eq. I|A2(1 . which are at least of third order in 
the marginal and relevant parameters, should not be in- 
cluded. However, there is a contribution proportional to 
vf arising from the eight-point vertex which must be in- 
cluded. To see this, let us explicitly write up to sec- 
ond order in the marginal and relevant parameters (and 
zeroth order in 77), 



r ; (6) (qi,q'2,q3;q3,q2,qi) = { 35 (i,2),3 (G ; (qi + q2 - q)(ttz + a/(|qi - q| + |q 2 - q|) + fy|qi + q2|)) q 

+3S(i', 2 '),3"%,2),3 (G;(q3 - q' 3 + q)(«z + ai(\q 3 - qg|) + 6;|q 3 + q|))^ + (q, <-> q-)} 

^ «fteH^- 8 )'(rj? t (e- ( ql 1 ... ) e- ( qJ,e-'<iie- < 4e- t q3,.. 1 e- t q 1 )) , (6.10) 
1 + n Jo \ ' q 

I 



where we used 

lfVi,...,qU4,...,qi) = f dt e-^-^ Vr Hs-mt 

Jo 

x t\%(e~%, . . . , e-*qi; e"^, . . . , e"*^) , (6.11) 



and from Eq. (| A3|) we have 

r| 8) (q'i, ■ ■ ■ ,q4;q4, ■ ■ ■ ,qi) = 
-2uf ji6S ( i !2 ,3),4<V,2%3'),4' »(q'i + q 2 + q3 - <u) 

+18<%,2),(3 ) 4)<%/ >2 '),(3',4') Xi(qi + q2 - q'i - q 2 )| 

+ terms at least cubic in v\ , a;, bi, and ui . (6-12) 
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The definition of the symmetrization operato rs £( 1,2), 3 
used in Eq. (|6.10|) and similar ones used in Eq. 1)6. 12[) can 
be found in Eqs. (|AlblA4alA4b|l of Appendix^ Keeping 
only the zeroth order in a momentum expansion, we have 

f[ 6) ({q, ; = 0}) - -12/W; - 15/W; - WoVih (6.13) 



dt e-^ D -^vU[mi-M + ISxi-tie-')] 



where, instead of Eqs. l|4.23l FO^ . we now have to cal- 
culate the zeroth and first order terms in a momentum 
expansion of Eq. (|6.15|l in powers of both the relevant and 
marginal parameters. We discuss the terms in Eq. (|6.15|) 
below, beginning with the contribution from the six-point 
vertex. 



1 + n J 

where 



0o = *i(0) = 



1 



2 ' 



(6.14) 



B. Improved truncation of the flow equation for 
the four-point vertex 

Including the terms which become marginal in D = 3 
we need to calculate 



r[ 4) (q' 1 ,q 2 ;q 2 ,qi) 



p(4,BCS)f„/ 

f,(4,ZS)/„/ 
1 I 



(qi,q 2 ;q2,qi) 



(qi,q 2 ;q2,qi) 



A(4,ZS')/ / 



(qi,q 2 ;q 2 ,qi) , (6.15) 



1. Contributions from f 



To calculate (f^ (q^, q' 2 , q; q, q 2 , qi))- , we write, as- 



(<■>) 



suming f ;f. = 0, 



fr ) (q' 1 , q2 ,q; q ,q 2 , qi ) = t dt e^'-* 

Jo 

x f ^(e-'q'i, e-%, e-% e'% e^q*, e^qx). (6.16) 



Using Eqs. (|6. 1616. 11)16. 121) . we can calculate the contri- 
bution of f[ 6 ^ to fl | . Including from the eight-point 
vertex only the q^ = contribution and expanding the 
other contributions up to linear order in the external mo- 
menta we have, 



r[ 6) (Aqi , Aq' 2 , Aq, Aq, Aq 2 , Aqi ; 



{l2«,C9b+Xj(A))+o,(l2 /„- S\,(A)(1 + A) h- n<I>,( A) + 2 y -,( A) 



+k (8/3 + 4x/(A)(l + A) + 4$,(A) + 2ft(A)) + A|q x + q 2 | (2ftuj + (40i + 2/3 2 )a, + 2/3 6 ; ) 
+A (|qi - qi| + |qi - q' 2 |) {4,0m + (4ft + 40x + 2/3 2 )a, + (4ft + 2/3 2 )6 ; ) } 

dt e-^ D ~ 8 Kl t [I6xi-t(0) + 18xi_ t (e-*)] , (6.17) 
I 



l + n 



where 



01 = Xi(0) = 



So 



(1 + r,) 2 (l + n) 3 ' 



02 = 



2S*i 



(l + n) 2 



(6.18) 
(6.19) 



and 5*0 and Si are defined in Eqs. 1)4.15)1 and 1)4.16)1 . We 
further introduced the functions 



(l + n) 5 



q + q'A 



(6.21) 



Since to zeroth order in A Eq. 1)6.17)1 determines the 
flow of vi via Eq. 1)6. 8fl . the integral on the right hand side 
of Eq. 1)6.16(1 gives simply vi to order A . We therefore 
split the angular average of the six-point vertex into two 
contributions, 



*/(A) = T -—(G i (q + q'A) |q + q'A|) , (6.20) 
l + n \ I q 



and 



<f[ 6) (qi ; q 2l q; q,q 2 ,qi)) a 



«i + (f I (w) (qi > qi > 4;4q3.qi)>a . ( 6 - 22 ) 
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where T\ contains the irrelevant parts of the six-point 
vertex. To calculate the angular average of the irrelevant 
part, we can neglect the rji dependence of the integral and 
also approximate the retarded dependence of the param- 



eters on I — t by an I dependence, since the retardation is 
exponentially damped. For large I we therefore can write 
the contributions of the irrelevant part to Eq. (j6.16(l as 
follows 



(f| 6l) (qi,q^q;q,q 2 , qi )). « J y{6 Ui fc(A) - O ) + a, (4*,(A)(1 + A) + 4$ ; (A) - 8/3 + fa (A) 

+h (2x/(A)(l + A) + 2$,(A) - 4/3 + \tpi{\) - &(0)]) + A|qi + q 2 | (0m + (20i + /3 2 )a ; + /5 60 
+A (| qi - qi| + | qi - q 2 |) (2/3i«, + {20o + 20 x + /3 2 ) a/ + (2ft + &)&,) } . (6.23) 



We now turn to the calculation of the BCS, ZS and ZS' 
diagrams in Eq. I|6.15|l . Keeping only the relevant and 
marginal terms in the expansion of the four-point vertex, 
see Eq. (|6.3ll . yields all contributions to second order in 
the marginal and relevant parameters arising from these 
diagrams and the irrelevant part of the four-point vertex 
does not contribute at this order. 



2. Contribution of the BCS, ZS and ZS 1 diagrams 



The contribution from the BCS-channcl to Eq. (|6.15(l 



is 



f,(4,BCS)/ / / 



I 



(qi,q 2 ;q 2 ,qi) ~ 
l 2 
~2 i Ul + 6 'lqi + q2|] Xidqi + qal) 

-a t [m + fy|qi +q 2 |] 

1 /6(1< |qi+q 2 -q| < e') 



l + n\ |qi + q 2 - q| 2 + rj 

x [|q - qil + |q - q 2 l + |q - q2| + |q - qil 

af /0(1< |qi+q 2 -q| < e') 



1 + rjX |qi + q 2 - q| 2 + n) 
:(|q'-q , il + |q'-q 2 |)(|q-q 2 | + |q-qi|))., (6.24) 



and the contribution from the zero-sound channel is 



rf- ZS) (q'i,q 2 ;q 2 ,qi) « 

-[ui +a ; |qi -qi|] 2 X;(|qi -q'll) 
-(a/ + hi) [m + a ; |qi - q' x |] 

1 /6(1< |qi-qi+q| < e l ) 



x ■ 
l + n 

x[|q-q'il 



|qi - qi 
|q + q 2 |- 



qr + n 
|q- q 2 | 



qil] 



2 , 2 , 1 /9(1< |qi-qi+q| < e 1 



<af + bf) 



1 + n \ |qi - qi + q|* + n 

[|q-q'il|q-q2| + |q + qi||q + q 2 |] 



1 /e(K|qa-qi + q|<e') 



1 + n \ |qi - qi + q|" + rj 

x[|q-q'il|q + q 2 | + |q + qi||q-q 2 | 



(6.25) 



The contribution from the other zero-sound channel 
(ZS' in Fig. ^) is obtained by replacing q' : <-* q 2 on the 
right-hand side of Eq. (|6.25l) . Expanding to linear order 
in the momenta, we arrive at 



r iin .i (qi,q 2 ;q 2 ,qi) 



-0o 



+ 2a, + 2uiai 



|qi + q 2 | 



1 .2 



uf + 2(/?i + 2 )af + {20x + 2 )uiai 



+0 o uibi + 2{3 aibi 



(6.26) 



Similarly, we obtain for the contribution from the zero- 
sound-channeL 



n(4,ZS) 



r 



(qi,q 2 ;q 2 ,qi) 



-0o [ui + at + bi] 



|qi - qil 



Pmf + (20o +0!+ 02)4 + (0! + 2 )bj 

+ {20o + 20! + 2 )u m + {20, + 2 )mbi 



+2{0 O + 0i + 2 )aik 



(6.27) 



The contribution from ZS' is again obtained by replac- 
ing qi <-> q 2 . 
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C. Truncated flow equation for the two-point 

vertex 

To calculate the flow of r; , we need to calculate T\ (0), 
see Eq. <|3.18[) . We write it as 



(l + n)f ; (2) (0) = (f ; (4) (0,q,q,0)). « e '- 2/ o dT "^o+ / (fte t - 2 J"«-**^(f l ( i ) t (0 > e- t ^e- t q J 0)). (6.28) 



with the second order expression of f , given by Eq. (|6.15|1 . Thus, we have 



f [ 4) (0, Aq, Aq, 0) « --Xi(X) [(uj + \k + ai f + (u, + h + a, A) 2 + (uj + a,[l + A]) 2 ] 
— ^(A)[oj(uj + 6; A + a;) + ai{ui + bi + a;A) + + ai[l + A])] 
-^z(A) [2a 2 + o 2 ] - («j + h + a t ) [fi a ui + (6, + aj)pi(A)] 

+ T ^<f[ 6) (0,Aq,q / ;q',Aq,0)> a( , 

1 + Tl X '1 



(6.29) 



with 



1 + n ; G Kq + q'A)|q + q'A|^ 



(6.30) 



Keeping only terms up to linear order in A, which define 
the flow of U[ , ai and 6/ , one finds 



A(2) rn x _ Ul + «i + 



(6.31) 



which is just the contribution of the relevant and 
marginal parts of ^ (Q, q, q, 0))^. Higher orders in A 
correspond to the contribution of the irrelevant part of 
the four-point vertex. With 

<f| 4) (0,q,q,0)). = Ui + ai+bi 

+ (f| 4i) (0,q,q,0)) a (6.32) 



D. Flow equations of marginal and relevant 
parameters 



It is now straightforward to write down the flow equa- 
tions for the marginal and relevant coupling parameters. 
The flow of vi is determined by Eq. JEBJ. With Eq. (|6~T3|) 
one finds 



d m = (6-2D-37]i)vi-3j3 Q vi(Aui + 5ai + 3b l ) 



l + riJ 



dX A 4D - 9 [8x/(0) + 9x/(A)] ,(6.34) 



we can write the contributions of the irrelevant parts (we 
here ignore the contributions of the irrelevant part of 
the six-point vertex beyond those which are implicitly 
contained in the renormalization of ui) as 



f>(0,q,q,0)>, « /^(r| 4 - BCS )(0,Aq;Aq,0) 
+lf ' zs) (0, Aq; Aq, 0) + f [ 4 < zs ' } (0, Aq; Aq, 0) 



_p(4,BCS) 



lni , (0,AfcA$,0)-rj&f ) (0,A&AQ,0) 

(6.33) 



-fgff >(0, Aq;Aq,0)\ , 



which is identical to the second and higher order in A 
contributions contained in Eq. 



where we approximated the retarded dependence on vi-t 
and xi-t in Eq. I|6.13[) by vi and \i (since the retardation 
is exponentially damped) and further took for the lower 
bound of the integral the limit I — > oo. Eq. (|6.34|l obvi- 
ously has au, =0 fixed point solution, which is stable if 
3/3 (4u* + 5a* + 36*) > 6 - 2D - 3r). 



Collecting 
from Eqs. 

r(4) 



all contributions 
(I6~2H l6~2^1 l6~27ll . 



to Eq. (|ETT5|) 

and decomposing 



T; (qi j qai q2) qi) into a momentum independent 
contribution and linear contributions according to 
Eq. EH, we find, using Eqs. j3~T31 IO the 

flow equations for the parameters characterizing the 
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four-point vertex, 



and 



dm = (4 - D - 2r]i)ui - -Porf - (3 ui(6ai + 4bi) 



-2(3 {2af + bf + 2a t bi) + 



vi 



1 + 77 l + n 



1 d\ 

T 



x |6 W; (xi (A) - O ) + ai (a X i (A) (1 + A) + 4$, (A) 
-80 o + [<MA) - 0/(0)]) + 6;(2»(A)(1 + A) 
+2$ ; (A)-4/?o + [^(A)-^(0)])} , (6.35) 

<9,a, = (3-Z?-2? 7; )a/-/3iu?-(2/3o + 2/3i+/3 2 )u ; a / 
-(20i + (3 2 )uibi - (2(3 + 0i + f3 2 )a 2 
-(ft 2(/% + ft +&)a,&, 



{2/3 lU; + (2/3 + 2/3 1 +/3 2 )a / 
+ (2/? 1 +/3 2 )fo ; } , 

01 2 



(6.36) 



{3-D- 277)6; - ^ttf - (2/?! + /3 2 ) Ui a, 

-0 o wi&z - 2(0i + 2 )a 2 - 2[3 aibi 
vi 



l + n 



{P 1 u l + (2(3 1 +/3 2 )a l +(3 b l } . (6.37) 



The flow equation for 77 is obtained from Eq. I|3.18[l and 
Eqs. (I6.28l6.33t . 



dm = (2 - m)n + 
l 



+ ai + 6; 
l + r ; 



;i hr -(f^(0,q;q,0)>. 



(6.38) 



Finally, we need to determine the flowing anomalous di- 
mension rji before we can analyse the fixed point of the 
RG equations. For now, we keep only the marginal and 
relevant part of the four-point vertex in Eq. I|3.21[) and 
arrive at 



m 



03 (aj + h) , 



where 



03 



d 



l + ri dq 2 



<|q' + q|)a< 



D — 1 

qi=a ~ 2D(l+n) 



(6.39) 
(6.40) 



1. Fixed point values in D — 3 

In D = 3, we have 

A + y/—Ti arctanh 



<MA) 



l+n+A 



A(i + n) 

A(2 + A) + rj In 



.41) 



(l+A) 2 +n 



2A(l + r ; ) 



3 + A 2 
3(1 +n) 2 ' 



(6.42) 
(6.43) 



0o 
0i 

02 
03 



u -i 


- n) 




- l 


2(1 


+ nf 




l 


2(1 


+ n) 2 




l 


3(1 


+ n) 



(6.44a) 
(6.44b) 
(6.44c) 
(6.44d) 



Setting the left-hand sides of the five flow equations 
(|6. 3416. 3516. 3616. 3 716. 38|) for the two relevant coupling 
parameters ri and ui and the three marginal coupling 
parameters a/, bi and vi equal to zero and employing 
Eq. (|6.39|) for the anomalous dimension, we obtain nu- 
merically the fixed point values 



r* « -0.0996 , « 0.122 
6* « 0.0339 , u* = , 



a* w 0.0371 , (a ak\ 
« 0.0263. (6 ' 45) 



The value of 77 is now much closer to the correct value 
77 w 0.038 and the inclusion of the marginal terms cer- 
tainly improves upon the analysis including only ui. 
However, the fact that = is an artifact of the ap- 
proximation which ignores all terms of third or higher 
order in the relevant and marginal parameters. A simple 
improvement can be obtained by including from the third 
order terms in the renormalization of vi the marginal con- 
tributions. In that case, only the flow of vi is modified 
and Eq. (|6.34|) becomes 

dm = {6-2D-3ni)vt-3p Vi(4;Ui + 5ai + 3bi) 



1 + n jo 

(ui +ai + b t ) 



(l + n) 3 



d\ \ iD - 9 [8 X i(0) + 9xiW] 
12( U/ + ai +b{f 



-9(uj + 2aif 



(6.46) 



The resulting fixed point values with the 3rd order terms 
are: 



r* w -0.134 , 
K « 0.0639 , 



: 0.127 , 
0.0917 . 



s 0.0693 
0.0513 . 



(6.47) 



The value for 77 is now slightly too large and a quick con- 
vergence is not obtained. However, one may expect from 
this result that a consistent treatment to third order, 
a very complex calculation, would indeed improve upon 
the second order result. For completeness, let us also 
mention the results of another possible approximation, 
where one keeps only the marginal and relevant terms of 
the vertices entering the f ^ expressions. Within such 
an approximation, the flow equations for a; and 6/ remain 
identical to Eqs. (|6.36I6.37|) whereas the flow equations 
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for ui, vi and n would simplify to 

dm = (6 - 2D - 3m)vi - 3/9b«j(4uj + 5a; + 36j) 



(it; + a ( + 6 ( ) 



(l + n) 3 



12{ui + ai + bif 
+9( Ul + 2 ai ) 2 



(6.48) 

<9 ; u ; = (4 -D- 2r)i)ui - ^/3 uf - /3 ui(6ai + 46/) 

(6.49) 

(6.50) 



-2/3 (2 a 2 + 6 2 + 2 a; 6 ; ) + 



i + n 



dm 



m + ai + h 



l + n 

The corresponding fixed point values are 

r* re -0.227 , u* « 0.178 , a* « 0.0838 , 
6* PS 0.0767 , v* ps 0.255 , 77 ps 0.0692 . 



(6.51) 



The value for rj is worse in this approximation compared 
to those which also include irrelevant parts of the vertices. 
Problems with a similar implementation of the sharp cut- 
off formulation were previously reported in Ref . [26j . Ig- 
noring the third order terms in the flow of vi, Eq. H6.48fl . 
leads to fixed point values almost identical to those listed 
in Eqs. ljQ5j) . 



2. Including irrelevant terms in the equation for n 



Eq. (|6.39|) for the flowing anomalous dimension in- 
cludes only the marginal parts of the four-point vertex. 
To include also irrelevant contributions we need to eval- 
uate the irrelevant parts of r[ 4 ^ , i. e. we must calculate 
both the irrelevant contribution from the six-point vertex 
and the BCS, ZS and ZS' contributions to 



1 + r, Jo 
d 

X ■— r 



dAA" 2+2 " 



f£)(Aq,Aq',Aq',Aq) 



(6.52) 



q 2 =0 



using Eqs. I|6. 1016. 1616. 2316. 2416. 2g|) . This is a somewhat 
lengthy calculation, and we refer to Appendix [B] for a 
table of the required integrals. The resulting values of 77 
are generally very small, in the approximation with all 
terms up to second order in the relevant and marginal 
parameters we obtain 77 ps 0.0127. Including the third 
order terms in the flow of vi we get the even smaller 
value 77 ps 0.00797. It is not completely clear why the re- 
sults for 77 become worse on including irrelevant terms in 
the calculation for r\i. A possible problem is that in the 
calculation of 77, Eq. I|6.52f) . irrelevant terms contribute 
which however do not also enter the renormalization of 
Ui, ai and 6/. Terms similar to the irrelevant ones con- 
tained in Eq. I|6.52|l would only enter the renormalization 
of other parameters at third order in the marginal and 
relevant parameters. This might suggest that one should 



use the same level of approximation for the vertex f [ 

in calculating 77 and for the flow equations for ui, 
ai and 6;. In that case, the first order expression for 77, 
Eq. (|6.39|) . should be used together with the second or- 
der expressions for the flows of ui, ai and 6;. As shown 
in the previous subsection, this indeed leads to better 
approximations of 77. 



VII. CALCULATION OF a(x) INCLUDING 
ONLY MARGINAL PARAMETERS 

We demonstrate here the importance of the irrelevant 
parts of the four-point vertex for a correct description 
of the large wave vector regime and show that a treat- 
ment that ignores irrelevant terms leads to an incorrect 
description of a(x) in the regime x 3> 1 and the wrong 
result AT C cx uqItiuq. Ignoring in Eq. I|6.3|) the irrele- 
vant contributions, one obtains from Eq. (|3.12|) the fol- 
lowing expression for the subtracted function f f 171 ' (q) = 
fp) ( ^_r« 



rj +0) (the superscript m indicates that only 



marginal terms remain) 



,(2m) 



(q) 



ai + bi 

i + n 

ai + h 

i + n 



(<|q' + qlW 



q - faq<l, 
V-3g+i f rg>l. 1 ' 

3q ^ 



Via Eq. (|4.49() we obtain the scaling function a(x) (we 
here restrict the discussion to D = 3), 



a x 



h ^ dl c -2»-u+firfT„ T a i + h 



3(l + r ; ) 



Ji a -\ux 3(1 + n) 

3x 



(7.2) 



This equation is valid for lnx < l c , for In a; > l c the 
integral bound l c — In a; is replaced by zero, i. e. the first 
integral vanishes. Since in D = 3 we have (o; + 6/ )/ (3 [1 + 
n]) ~ Vi ( see Eq. H6.39JI ). o-(x) is determined by the flow 
of rji alone in this approximation. Let us examine first 
the case x ^> 1 . The first integral then extends only over 
small I <^ l c . In that case, all marginal parameters are 
small and we may ignore the contribution of the integral 
of i] T in the exponent and approximate the flow of the 
marginal parameters as a;, 6/ tx e 2 ^ l ~ lc ^u1 since 9;a; w 
29(6; « uf oc e 2 ^ l ~ lc ^u 2 for I <g; l c . The contribution 
of the first integral in Eq. I|7.2|l therefore vanishes like 
x 2 (l c — lnx)ul as x approaches the UV cutoff, x — > e lc 
and makes no contribution for x > e' c . To investigate 
the contribution of the second integral, we split it into 
two parts, 



-In s 



(7.3) 
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where, to estimate the contribution of the regime I < l c , 
we use the same approximation as for the first integral 
in Eq. 1)7. 2J1 . This yields a contribution oc (x — 1 — In x — 
(1 — x)/3x). To estimate the contribution of the inte- 
gral for I > l c we may replace the parameters a;, bi, ri 
and rji by their fixed point values and obtain a contri- 
bution oc rix/(l — r/) — r\j(2 — rf) + rj/(3x(3 — ??)). We 
have employed Eq. (|(j.39|) to arrive at this result. Hence, 
the dominant behavior of cr{x) for large x is linear in x. 
This large x behavior prohibits the removal of the UV 
cutoff in the integral determining the temperature shift, 
Eq. 15. and one thus would incorrectly predict from 
Eq. (|5.1I) a term for the T c shift behaving like uq In uq in 
D = 3. Thus, inclusion of the irrelevant terms is essential 
for determining correctly the shift of the critical temper- 
ature. The absence of irrelevant terms in the Wilson RG 
analysis carried out by Bijlsma and Stoof [t| seems to 
be responsible for the uq In uq behavior they find for the 
T c -shift. 

On the other hand, an approach which explicitly in- 
cludes marginal terms but ignores irrelevant ones is well 
suited to describe the scaling regime i->0. In this case, 
the dominant contributions to Eq. (|7.2|l come from I > l c 
where we replace the coupling constants by their fixed 
point values and obtain a(x) ~ A 3 x 2 ~ n with 

* = i + A-2^^" 1 + T*- (7 ' 4) 

1 — rj 2 — rj o — -q 6 

This expression is in relatively good agreement with our 
result from Sec. IIV D 21 Eq. I|4.57[l . where we obtained 
A 3 ps 1.17 for r\ ps 0.104, whereas Eq. (T7T^f) would predict 
A 3 ps 1.22 for the same 77. Furthermore, for the generally 
accepted value rj ps 0.038 Eq. ifTIJ) predicts A 3 ps 1.07, 
which compares relatively well with a recent Monte-Carlo 
result A] A 3 ps 1.04. 

VIII. CONCLUSION 

In this work we have shown how the functional RG 
can be employed to calculate the complete scaling func- 
tion of the zero frequency self-energy of weakly inter- 
acting bosons in 3 < D < 4. The scaling function de- 
scribes the cross-over from the small wave vector regime 
with anomalous scaling to the large wave vector regime 
where logarithmic divergences appear in D = 3. A simple 
truncation of the flow equations at the four-point vertex 
which ignores the six-point vertex and does not treat the 
marginal parts of the four-point vertex consistently leads 
nonetheless to a very accurate description of both the 
cross-over regime and the large wave vector regime while 
giving satisfying results for the anomalous scaling regime 
even in D = 3. We have used this scaling function to cal- 
culate the interaction induced shift of the critical temper- 
ature and obtained AT C /T C = 1.23 a?! 1 / 3 , a result which 
compares very well with those obtained within the varia- 
tional perturbation theory (l6lll7| and Monte-Carlo sim- 
ulations [lil n~5| . The technique and truncation scheme 



employed here might also prove useful to obtain energy 
and/or momentum-dependent scaling functions in other 
critical systems. A similar truncation has already been 
applied to calculate the self-energy in the vicinity of the 
Luttinger liquid fixed point of fermions in one dimension 

m 

We have further investigated the fixed point structure 
using several approximation schemes which include the 
marginal terms associated with the four-point and six- 
point vertices in D = 3. These schemes generally lead 
to an improved anomalous dimension, however, we did 
not get a quick convergence of 77 to the accepted value. 
Nonetheless, the best value for rj which we obtain, rj ps 
0.0513 is rather close to the one obtained from the first 
order average action approximation, -q ps 0.049 j2^ and it 
is feasible that an improved treatment, e. g. a consistent 
treatment to third order in the marginal and relevant 
parameters, would indeed produce quite accurate results. 
Unfortunately, such a calculation is rather lengthy and 
seems not to be an efficient way of calculating accurate 
values of critical exponents. 

Finally, we have analysed the two-point scaling func- 
tion within an approach which includes marginal and ig- 
nores irrelevant terms of the four-point function. It was 
shown, that if one ignores the irrelevant terms contained 
in the higher-order momentum dependence of the four- 
point vertex, one obtains an incorrect UV behavior of 
the scaling function and hence the wrong functional de- 
pendence AT C /T C oc an 1 ' 3 lnfan 1 ' 3 ) of the critical tem- 
perature shift in D = 3. However, the approach worked 
well in the critical regime where we used it to calculate 
the prefactor of the anomalous scaling term, which was 
shown to be in good agreement with numerical results 
■ While we have not attempted to calculate the scaling 
function including both marginal and irrelevant terms, it 
seems certainly feasible to do so within an approximation 
which ignores the six-point vertex (since v* = in a cal- 
culation to second order in ui, a;, bi and vi, this would be 
consistent). The fixed point value of ij and hence the de- 
scription of the small wave vector regime improves within 
such an approach. However, we do not expect that the 
inclusion of marginal terms would have much effect on 
the crossover regime x ps 1 , since this regime seems to be 
well described already in the simpler truncation used in 
Sec. II VI as is evident from the rather accurate value for 
the T c -shift. 



APPENDIX A: FLOW EQUATIONS FOR THE 
SIX- AND EIGHT-POINT VERTEX 

Below, we discuss the flow equations of the six- and 
eight-point vertex for the unrescaled vertex functions. 
The inhomogeneous part of the flow equations for the 
classical rescaled flow equations follow via the replace- 
ment rule G A (K) -» -G ; (q), G A (K) -» -G,(q), 
p(2n) _^ p( 2 ™)^ _i j anc j multiplying the resulting 
expression by an overall minus sign. 



21 



1. Six-point vertex 

For the calculation of the critical exponents in D = 
3 we need the six-point vertex, see Sec. IVII We first 
define the symmetrization operators £1,2,3 and £1,(2.3) a $ 
follows, 

5i, 2 , 3 /(l,2,3) = i[/(l,2,3) + /(2,3,l) + /(3,l,2) 
+/(3,2,1) + /(2,1,3) + /(1,3,2)] , (Ala) 



Given a function /(l, 2, 3) that is already symmetric with 
respect to the pair (2,3), the function 5w 2 ,3)/(l> 2, 3) is 
a totally symmetric function. The flow equation of the 
six-point vertex is given by |13| (the equation is shown 
graphically in Fig. [SJ 



5i, (2 , 3) /(l,2,3) - 5 (2 ,3),i/(l,2,3) 

= i [/(l, 2, 3) + /(2, 1, 3) + /(3, 2, l)IAlb) 



J 



+3 / { £3,(2,1) fGAWG A (if')ri 6) (^UU3; K z , K',K)Tf{K,K'; K 2 ,K 1 )} 

+<5(i', 2 '),3' [GA(^)G A (if')ri 4) (#i\ ^2; ^)ri 6) (i^, K', K' 3 ; K 3 ,K 2 ,K 1 )} } 
+9 ^(^^.(^[[Ga^ 

+9 y , 2 '), 3 '<S3,(2,i) [[G a (K)Ga(K')G a (K") + Ga(K)G a (K')Ga{K") + G A (K)G A {K')G A {K")] 

xrf(K[,K 2 ; K, K')r^(K' 3 ,K'; K" , K 3 )vf{K", K; K 2 ,K : )\ K "^ +K ^ K 
+36 / 5i/,a/,s/5i A 3 [G A (^)G A (^')GA(^")r A 4) (^, i^'; AT, Ki)^^, if"; K', Jfajr^^, K; K" , K 3 )] * =Ki ~* : 



VK 2 -K 



K"=K' 3 -K 3 +K 
K 

(A2) 



2. Eight-point vertex 

In Eq. l|tj.lU[l we also need the following terms from the 
inhomogeneity of the eight-point vertex, 



d A r^(K' u K' 2 , K' 3 , K' A ; K 4 , K 3 ,K 2 , K x ) 

= 16 / «S (1 , 2 , 3) ,45 (1 ,, 2 ,,3').4'[ G A(^) G A(^')rf (K' u K , a ,K' 3 ;K A ,K',lC)T^(K,K',K , 4 ;K !l ,K a ,K 1 )\ 

Jk L 1 k'=k' 1 +k 2 +k! ) -k 4 -k 

+9 J^S (1 , t2l)t(3 , At) S {lt2)i{3t4) [[GA(K)GA(K') + Ga(K)G a (K')] 

xTfiKlK^K'-K^^K^vfiK^K^K-K',^,^)] +..., (A3) 

I 

;h 5 (1 , 2)i(3 ,4)/(l,2,3,4) = i[/(l,2,3,4) + /(l,3,2,4) 

5 ( i,2,3), 4 /(l,2,3,4) = i[/(l,2,3,4) + /(4,2,3,l) +/(1, 4, 2, 3)] . (A4b) 

+/(1,4,3,2) + /(1,2,4,3)] , (A4a) 
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A: 



36 




FIG. 6: Diagrammatic representation of the flow equation 
for the six-point vertex, see Eq. 1A2I . 



APPENDIX B: SOME INTEGRALS ENTERING 
THE CONTRIBUTION FROM IRRELEVANT 
TERMS TO r\ 

We here calculate the irrelevant contributions to the 
flow of 77, using Eq. llfT52|) of Sec. IVID2I Let us define 



//(qi,q2) = 



.9;(qi,q2) 



I /9(1< |qi+q 2 -q| < e l ) 



l + n\ |qi+q 2 -q| 2 
x [|q- qi| + |q- qal 



n 



(Bl) 



1 /8(1 < |qi +q 2 -q| < 



1 + n \ |qi + q 2 - q 



ft|2 



[|q-qil 



q- qa 



(B2) 



+ t^\t + 

4(f )(> 



FIG. 7: Diagrammatic representation of the terms of the 
flow equation for the eight-point vertex which are of order vf , 
see Eq. <A3t . The omitted terms are at least cubic in the 
parameters ai,bi,vi and it;. 



J 



dq 2 



rf(Aq,Aq',A'q";A'q",Aq',Aq) 



q',q" 



q 2 =0 



hi(q.i,q. 2 ) = 



2 /6(1< |qi+q 2 -q| < e l ) 



l + n 

x|q-qiliq-q2 



qi + q 2 - q| 2 + n 



(B3) 



To calculate the contributions of the six-point vertex 
to Eq. I|6.52|) . we need to evaluate (we restrict the results 
to D = 3) 



-«,{3u, [x'{(\, A) + xi'(\ A')] + 2(6 ; + 2o,)[/,"(A, A) + //'(A, A') + d,(A, A')] + (b t + oi)^'(A)} , (B4) 



with 



//'(A, A') = A (/;(Aq,A'q 



x;'(A,A') 
di(A,A') 



= X/(0)A(|q' + Aq| 



_ -A 2 (l + A') 2 + A 2 (3 + 2A')r ; 
, 2=0 ~ 12A'(l + r ; )[(l + A') 2 +r ; ] 2 
A 2 



5 s=o 3(1 + r ; ) 2 



_2 9_ 

l + n <9g 2 
_L d_ 

1 + ndq* 



Gi (Aq + A'q" + q') ) 

/ q',q'' 



Gi (Aq + A'q" + q')|Aq + A'q : 



A 2 [ n -(i + A') 



A21 



q2=0 6(i + n)A'[(i + A') 2 + n]2 ' 

_ A 2 [(l + A') 2 (2 + A') + (2 + 3A')n] 

q 2 =0 



12(l + r ; )A'[(l + A') 2 + r ; ] 



(B5) 
(B6) 
(B7) 
(B8) 



To calculate the contributions of the BCS, ZS and ZS' channel to Eq. (|6.52() . further averages are required: 



q^2 (91 (Aq,Aq')) q , 



A(l + A) 2 (2 + 3A) + (-4 + (-2 + A)A)r ; + 2((1 + A) 2 + n f In ( (1+ 1 + ) ? +n ) 



q 2 =0 



i2(i + n )[(i + A) 2 + nr 



, (B9) 
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dq 2 



y,(Aq, Aq')|Aq + Aq' 



q' 



A [8(1 + n) + A [(3 + A) 2 (3 + A(2 + A)) + 2(10 + A(6 + A))n + rj 



24(l + r,)[(l + A) 2 +r ; ] 



2(l-A 2 + n ) [(l + A) 2 +r ; ] 2 arctanh(^^) +(l + A 2 +n) [(l + A) 2 +r ; ] 



24(l + r ; )[(l + A) 2 + n ]'y^ 



(BIO) 



^2 ( Af i> A q" 



q 2 =0 



[(A - 1) (1 + A) 2 (A + A 2 + 2r,) + Ar 2 ] ^ + (1 - A 2 - r,) [(1 + A) 2 + r,] 2 arctanh 



6(l + rO[(l + A) 2 +r ; ] 



. (BH) 



xKO)^(|q" + Aq| | q" + Aq' | ) . , . „ 
d 



= A 2 (3 + A 2 ) 

g2 =0 ~ 9(1 + r; )2 



(B12) 



Xi(|M' + Aq|) |Aq' + Aq|< 



q' 



q 2 =0 



A 2 (l + A) 2 (4 + 3A) + A 2 (4 + 5A)n + A [(1 + A) 2 + n f In ( (1 + A + ) ', +? " ) 



6(l + r ; )[(l + A) 2 + r^ 



(B13) 



Using these expressions to calculate f $ from Eqs. ^TM3M2M2MM 

one can determine r\ via Eq. 1)6.520 . One 

can check that to linear order in A one has 



^(f j 4) (Aq, Aq', Aq', Aq)> \ q 2 =0 =s A(9 ;a; + %)/3 + 2A(a ; + 6,)»»/3 , 



(B14) 



which, after performing the A- integration, reproduces at the ^fixed point the result Eg. (|6.39|). 
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